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Abstract 

We devise a lucid and straightforward way for explicitly constructing quasi-periodic 
wave solutions (also called multi-periodic wave solutions) of supersymmetric equations in 



. superspace M^'^ over two-dimensional Grassmann algebra Gi{a). Once a nonlinear equa- 

\ tion is written in a bilinear form, its quasi-periodic wave solutions can be directly obtained 

by using a formula. Moreover, properties of these solutions are investigated in detail by 
analyzing their structures, plots and asymptotic behaviors. The relations between the 
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. quasi-periodic wave solutions and soliton solutions are rigorously established. It is shown 



that the soliton solutions can be obtained only as limiting cases of the quasi-periodic 
' 2 1 

' wave solutions under small amplitude limits in superspace M^' . We find that, in contrast 

to the purely bosonic case, there is an interesting infiuencing band occurred among the 

^ , quasi-periodic waves under the presence of the Grassmann variable. The quasi-periodic 

H ' 

I waves are symmetric about the band but collapse along with the band. Furthermore, 

the amplitudes of the quasi-periodic waves increase as the waves move away from the 
band. The efficiency of our proposed method can be demonstrated on a class variety of 
supersymmetric equations such as those considered in this paper, M = 1 supersymmetric 
KdV, Sawada-Kotera-Ramani and Ito's equations, as well as A/" = 2 supersymmetric KdV 
equation. 
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1. Introduction 

The algebro-geometric solutions or finite gap solutions of nonlinear equations 
were originally obtained on the KdV equation based on inverse spectral theory and 
algebro-geometric method developed by pioneers such as Novikov, Dubrovin, Mck- 
ean, Lax, Its and Matveev et al. [l]-[5] in the late 1970s. In fact, such a solution is an 
expression written in terms of the Riemann theta functions. Hence it is also called 
a quasi-periodic solution due to the quasi-periodicity of the theta functions. By 
now this theory has been extended to a large class of nonlinear integrable equations 
including sine-Gordon equation, Camassa-Holm equation, Thirring model equation, 
Kadomtsev-Petviashvili equation, Ablowitz-Ladik lattice and Toda lattice [6]-|16]. 

The quasi-periodic solutions have important applications in physics. For in- 
stance, they can describe the nonlinear interaction of several modes. All the main 
physical characteristics of the quasi-periodic solutions (wave numbers, phase veloci- 
ties, amplitudes of the interacting modes) are defined by a compact Riemann surface. 
There are numerous applications of the finite-gap integration theory in condensed 
matter physics, state physics and fiuid mechanics. For example, in peierls state, 
phonon produce a finite-gap potential for electrons, and the peierls state is a lattice 
of solutions at low densities of electrons [6]. A most famous mechanical system, 
the Kowalewski top, was the focus of interest in the 19th century. The equation of 
motion of the top can be solved through finite-gap theory |6J. A problem of fun- 
damental interest in fiuid mechanics is to provide an accurate description of waves 
on a water surface. The Kadomtsev-Petviashvili equation is known to describe the 
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evolution of waves in shallow water and admits a large family of quasi-periodic so- 
lutions. Each solution has independent phases. Experiments demonstrate the 
existence of genuinely two-dimensional shallow water waves that are full periodic in 
two spatial directions and time. The comparisons with experiments showed that the 
two-periodic wave solutions of the KP equation describe shallow water waves with 
much accuracy [T71 118] . 

The algebro-geometric theory, however, needs Lax pairs and involves compli- 
cated calculus on the Riemann surfaces. It is rather difficult to directly determine 
the characteristic parameters of waves such as frequencies and phase shifts for a 
function of given wave-numbers and amplitudes. On the other hand, the bilinear 
derivative method developed by Hirota is a powerful approach for constructing ex- 
act solution of nonlinear equations. Once a nonlinear equation is written in bilinear 
forms by a dependent variable transformation, then multi-soliton solutions are usu- 
ally obtained |19]-|25]. It was based on Hirota forms that Nakamura proposed a 
convenient way to construct a kind of quasi-periodic solutions of nonlinear equa- 
tions [261 EH EH], where the periodic wave solutions of the KdV equation and the 
Boussinesq equation were obtained. Such a method indeed exhibits some advan- 
tages over algebro-gometric methods. For example, it does not need any Lax pairs 
and Riemann surface for the considered equation, allows the explicit construction 
of multi-periodic wave solutions, only relies on the existence of the Hirota's bi- 
linear form, as well as all parameters appearing in Riemann matrix are arbitrary. 
Recently, further development was made to investigate the discrete Toda lattice, 
(2-|-l)-dimensional Kadomtsev-Petviashvili equation and Bogoyavlenskii's breaking 
soliton equation [29]- [M]. Indeed there are some differences between quasi-periodic 
solutions and algebro-geometric solutions. A quasi-periodic solution needs not be 
an algebro-geometric one. Sometimes a quasi-periodic solution may not correspond 
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to any Riemann surface and is generically associated with infinite bands, not just 
finitely-many, for instance with a Riemann surface of infinite genus. 

The concept of supersymmetry was originally introduced and developed for ap- 
plications in elementary particle physics thirty years ago [35]-[37]. It is found that 
supersymmetry can be applied to a variety of problems such as relativistic, non- 
relativistic physics and nuclear physics. In recent years, supersymmetry has been 
a subject of considerable interest both in physics and mathematics. The mathe- 
matical formulation of the supersymmetry is based on the introduction of Grass- 
mann variables along with the standard ones [38] • In a such way, a number of 
well known mathematical physical equations have been generalized into the super- 
symmetric analogues, such as supersymmetric versions of sine-Gordon, KdV, KP 
hierarchy, Boussinesq, MKdV etc. [39]-|l9]. It has been shown that these super- 
symmetric integrable systems possess bi-Hamiltonian structure, Painleve property, 
infinite many symmetries, Darboux transformation, Backlund transformation, bi- 
linear form and multi-soliton solutions. The systematic bilinear transcription of 
supersymmetric equations was introduced by Carstea [ISj HS]. This required an 
extension of the Hirota's bilinear operator to supersymmetric case. Despite this 
bilinearization of supersymmetric equations, the standard construction did not lead 
to malti-soliton solutions. In recent years, Carsta, Liu, Ghosh et al. have done 
much on the construction of soliton solutions of supersymmetric equations |l2]-[39] 
. However, the quasi-periodic solutions of the supersymmetric systems, which can 
be considered as a generalization of the soliton solutions, are still not available (both 
by algebro-geometric method and by bilinear methods or others) to the knowledge 
of the author. 

The motivation of this paper is to show how the quasi-periodic wave solutions 
of nonlinear supersymmetric equations can be constructed with Hirota's bilinear 
method in superspace. To achieve this aim, we devise a Riemannn theta function 
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formula, which actually provides us a lucid and straightforward way for applying in a 
class of nonlinear supersymmetric equations. Once a nonlinear equation is written in 
bilinear forms, then the quasi-periodic wave solutions of the nonlinear equation can 
be obtained directly by using the formula. This method considerably improves the 
key steps of the existing methods, where repetitive recursion and computation must 
be preformed for each equation [29]-[3lj. As illustrative example, we shall construct 
quasi-periodic wave solutions to the M = \ supersymmetric Sawada-Kotera-Ramani 
equation and M = 2 supersymmetric KdV equation. 

The organization of this paper is as follows. In section 2, we briefly give some 
properties on superspace and super-Hirota bilinear operators. In section 3, we intro- 
duce a super Riemann theta function and discuss its quasi-periodicity. In particular, 
we provide a key formula for constructing periodic wave solutions of supersymmetric 
equations. As applications of our method, in section 4 and section 5, we construct 
one- and two-periodic wave solutions to the M = 1 supersymmetric Sawada-Kotera- 
Ramani equation and Af = 2 supersymmetric KdV equation, respectively. The 
propagation of the quasi-periodic waves are displayed with help of software Mathe- 
matica. In addition, we further present a simple and effective limiting procedure to 
analyze asymptotic behavior of the periodic wave solutions. It is rigorously shown 
that the quasi-periodic wave solutions tend to the soliton solutions under small 
amplitude limits. At last, we briefly discuss the conditions on the construction of 
multi-periodic wave solutions of supersymmetric equations in section 6. 

2. Super space and super-Hirota bilinear form 

To fix the notations and make our presentation self-contained, we briefiy re- 
call some properties about superanalysis and super-Hirota bilinear operators. The 
details about superanalysis refer, for instance, to Vladimirov's work [SOj IST]. 
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A linear space A is called Z2-graded if it represented as a direct sum of two 
subspaces 

A = Ao © A, 

where elements of the spaces Aq and Ai are homogeneous. We assume that Aq is a 
subspace consisting of even elements and Ai is a subspace consisting of odd elements. 
For the element / G A we denote by /o and /i its even and odd components. A 
parity function is introduced on the A, namely, 

I I _ / 0, if / e Ao, 
1^1 -\ 1, if feA^. 

We introduce an annihilator of the set of odd elements by setting 

^Ai = {A e A : AAi = 0}. 

A superalgebra is a Z2-graded space A = Aq © A in which, besides usual oper- 
ations of addition and multiplication by numbers, a product of elements is defined 
with the usual distribution law: 

a{ab + /3c) — aab + (3ac, {ab + (3c)a = aba + (3ca, 

where a,b,c G A and a, 13 G C. Moreover, a structure on A is introduced of an 
associative algebra with a unite e and even multiplication i.e., the product of two 
even and two odd elements is an even element and the product of an even element 
by an odd one is an odd element: \ab\ — \a\ + \b\ mod (2). 

A commutative superalgebra with unit e = 1 is called a finite-dimensional Grass- 
mann algebra if it contains a system of anticommuting generators aj,j — 1, • • • ,n 
with the property: ajak + CkCTj — 0, j,k — 1,2, ■■■ ,n, in particular, crj = 0. The 
Grassmann algebra will be denote by G„ = Gn{(Ji, ■ ■ ■ , cr„). 

The monomials {eo, = aj^ ■ ■ ■ ctj^}, j — {ji < ■ ■ ■ < jn) form a basis in the 
Grassmann algebra G„, dimG„ = 2". Then it follows that any element of G„ is a 



linear combination of monomials aj^ • ■ ■ Uj^, ji < ■ • ■ < jk, that is, 

/ = /o + ^ X] fh-3k^h ■ ■ ■ ^3k ' 
fc>Oii<-<ifc 

where the coefficients fj^-j^. € C. 

Definition 1. Let A = Aq ® A be a commutative Banach superalgebra, then 
the Banach space 

is called a superspace of dimension (m, n) over A. In particular, if Aq = C and 
Ai = 0, then M™''' = C". 

A function f{x) : R^'" — >■ A is said to be superdifferentiable at the point 
X G R^'", if there exist elements Fj{x) in A, j = 1, • • • , m + n, such that 

m+n 

f{x + h) = fix) + J2 (^j(^)' hj) + o{x, h), 

3=1 

where x — (xi, ■ ■ ■ , Xm, Xm+i, ■ ■ ■ , Xn) with components Xj,j — 1, - ■ ■ ,m being even 
variable and Xm+j = Oj,j = 1, • • • ,n being Grassmann odd ones. The vector h = 
{hi, ■■■ ,hm, /im+1, • • • , h^+n) with (/ii, ■■■ ,hm) eK^ and (/i^+i, • • • , h,ri+n) e h-l- 
Moreover, 

II oix, h) II 
lim " , ; ,/ " — > 0. 

\\h\\^o II h II 

The Fj{x) are called the super partial derivative of / with respect to Xj at the point 
X and are denoted, respectively, by 

df{x) 



dxj 



The derivatives with respect to even variables Xj, j = 1, 2, • • - n are uniquely 



defined. While the derivatives to odd variables 9j = Xj^n, j = 1,2, •••m 



"J 



are not uniquely defined, but with an accuracy to within an addition constant 
C(7i • • • (7„, c e C from an annihilator ^Gn of finite-dimensional Grassmann algebra 



The super derivative also satisfies Leibniz formula 



= + [-imixf-^, i = 1, . . . , ™ + (2.1) 

Denote by P(A", A) the set of polynomials defined on with value in A. We say 
that a super integral is a map / : 7'(A^, A) — > A satisfying the following condition 
is an super integral about Grassmann variable 

(1) A linearity: /(/x/ + ug) = + ul{g), ^,ueA, f,ge V{A^, A); 

(2) translation invariance: /(/j) — /(/), where = f{6 + ^) for all ^ G A", 

/eP(A^A). 

We denote /(^^) = where s belongs to the set of multiindices Nn = {e — 
(£i, • • • ,en),£j = 0, 1, 0^ = 61^ ■ ■■ 61^ ^ 0}. In the case when = 0, £ e -/V„, \e\ < 
n — n — 1, such kind of integral has the form 

/(/) = J(/)/(l, •••,!), 

where 

Af)- 



Since the derivative is defined with an accurcy to with an additive constant form the 
annihilator -^L^, — {9i - ■ ■ 9^ e A"}, it follows that J : V ^ A/-'-L„ is single- 
valued mapping. This mapping also satisfies the conditions 1 and 2, and therefore 
we shall call it an integral and denote 

J{f) = J mde = Je,--- e^de, ■ ■ ■ do^, 

which has properties: 

^1 ■ ■ ■ ^nddi ■ ■ ■ d9, 

df 

^dOi-'-den^O, j ^l,--- ,n. (2.2) 
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In this paper, we consider functions with two ordinary even variables x, t and 
a Grassmann odd variable 6. The associated space M^'^ = Aq x Ai (we may take 
Aq = M or C) is a superspace over Grassmann algebra Gi(cr) = Gi^o © whose 
elements have the form 

f = fo + ficy. 

where e = 1 is a unit, a is anticommuting generator. The monomials {l,o"} form 
a basis of the G'i(cr), dimGi (cr) = 2. Therefore, any yU G Gi^i have the form = 
/3o", /3 G C. Under traveling wave frame in space ffi^'^, the phase variable should 
have the form 

^ = ax + ujt + [30a. 

For the functions f{x,t,6),g{x,t,6) : R^'^ — > A, the Hirota bilinear differential 
operators and Df about ordinary variables x, t are defined by 

D^D^fix, t, 9) ■ g{x, t, 9) = {d, - d^>r{dt - deTfix, t, 9)g{x', t', ^) U'=.,t'=t. 

The super-Hirota bilinear operator is defined as |12] 

D^-V(a:,t,^)2)^^7(x,t,^), 

where the differential operator D = de + 9dx is the super derivative, and the super 
binomial coefficients are defined by 

([J/|),if {N,j)y^ {0,1) mod 2, 

0, otherwise 
[k] is the integer part of the real number k {[k] < k < [k] + 1). 

We point out here that throughout this paper the natural number (which will 
denote powers, the number of phase variables, number of terms etc.) is different 
form A/" which is related to supersymmetry or superspace. 

Proposition 1. Suppose that functions f{x,t,9),g{x,t,9) : M^'^ — t- A, then Hi- 
rota bilinear operators D^, Dt and super-Hirota bilinear operator have properties 



s'!f(x.t,e)-g(x,t.e) 



A — n 



-1) 



J 
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m 

S'ff ■ 9 = D^f . g, 

D^D'^e^' ■ = {ai - a2)"'(wi - 002^6^'+^^ 



where = ajX + Ujt + f3j9a + Sj, aj, uj, aj, 6j G Aq are parameters, j = 1, 2. In fact, 
the third formula above is defined with an accuracy to within an addition constant 
of the ccr G"*- Ai. More generally, we have 

F{S„ D,, Dt)e^' ■ e^^ = F{a{Pi - P2) + - «2), «i - «2, c^i - c^2)e«^+^^ 

(2.3) 

where F{St, D^, Dt) is a polynomial about operators St,Dx and Df. This proper- 
ties are useful in deriving Hirota's bilinear form and constructing the quasi-periodic 
wave solutions of the supersymmetric equations. 

3. Super Riemann theta function and addition formulae 



In the following, we introduce a multi-dimensional super Riemann theta function 

2 1 

on superspace M^' and discuss its quasi-periodicity, which plays a central role in the 
construction of quasi-periodic solutions of supersymmetric equations. The multi- 
dimensional Riemann theta function reads 

£, s|r) = ^ exp{27ri(^ + £, n + s) — 7r(T(n + s), n + s)}. (3.1) 

Here the integer value vector n = (ni, ■ ■ ■ ^nj^Y' ^ 1 complex parameter vectors 
s = (si, ■ ■ ■ ,sj\f)'^,£ = (^1, ■ ■ ■ ,£n)'^ ^ C^. The complex phase variables ^ = 
(Ci; ■ • ■ ; ^nV, = OLjX + ujjt + (5j9a + 5j, Uj, Uj, l3j, 6j G Ao, j = 1, 2, ■ ■ ■ , A^, where 
X, t are ordinary variables and 6 is Grassmann variable. Moreover, for two vectors 
/ = (/i, ■ ■ ■ , /at)^ and g = {gi, ■ ■ ■ , gN)'^, their inner product is defined by 

(/, 9) = fl9l + /2^2 + ■ ■ ■ + fN9N- 
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The T = {Tij) is a positive definite and real- valued symmetric N x N matrix, which 
is independent of 6 and a in superspace M^'^. The entries Tij of the period matrix r 
can be considered as free parameters of the theta function (3.1). 

In this paper, we take the r to be pure imaginary matrix to make the theta 
function (3.1) real-valued. In the definition of the theta function (3.1), for the case 
s = £ = 0, hereafter we use t) = 0,0\t) for simplicity. Moreover, we 
have £, 0|r) = + e, r). It is obvious that the Riemann theta function (3.1) 

2 1 

converges absolutely and superdifferentiable on superspace M^' . 

Remark 1. The period matrix r here is different form algebro-geometric the- 
ory discussed in P-[33], where it is usually constructed via a compact Riemann 
surface F of genus N E N. We take two sets of regular cycle paths: ai, 02, ■ ■ ■ , cat; 
61, 62, ■ ■ ■ , &Af on r in such a way that the intersection numbers of cycles satisfies 

flfc o aj =bkO bj = 0, Ofc o bj = Skj, k,j = l,---,N. 

We choose the normalized holomorphic differentials ujj,j = 1, ■ ■ ■ , iV on F and let 

djk = ^j, bjk = ujj, 

•J J bk 

then N X N matrices A = (ajk) and B = (bjk) are invertible. Define matrices C 
and T by 

C = {c,k) = A-\ T = {r,k) = A 'B. 

It is can be shown that the matrix r is symmetric and has positive definite imaginary 
part. However, we see that the entries in such a matrix r are not free and difficult 
to be explicitly given. □ 

Definition 2. A function g{x,t) on x C is said to be quasi-periodic in t 
with fundamental periods Ti, ■ ■ ■ , G C if Ti, ■ ■ ■ , are linearly dependent over 
Z and there exists a function G{x,t) G x C'^ such that 

Gix, yi, ■ ■ ■ , Vj+Tj, ■ ■ ■ , 2/fc) = G{x, 2/1, ■ ■ ■ , ■ ■ ■ , 2/fe), for all G C, j = 1, ■ ■ ■ , k. 
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G{x, t, - ■ ■ 



t. 



,t) = g{x,t). 



In particular, g{x,t) becomes periodic with T if and only if Tj = mjT. □ 

Let's first see the periodicity of the theta function r). 

Proposition 2. [32] Let ej be the j— th column of N x N identity matrix 1^; 
Tj be the j— th column of r, and tjj the (j, j)-entry of r. Then the theta function 
r) has the periodic properties 



The theta function r) which satisfies the condition (4.4) is called a mul- 
tiplicative function. We regard the vectors {sj, j = 1, ■ ■ ■ , A^} and {irj, j = 
1, • • • , A^} as periods of the theta function ^{$,, r) with multipliers 1 and exp(— 27ri^j + ttTj-j), 
respectively. Here, only the first vectors are actually periods of the theta function 



Proposition 3. Let Sj and Tj be defined as above proposition 2. The mero- 



■i?(C + Bj + iTj, r) = exp{-2TTi^j + nTjj){}{^, r). 



(3.2) 




morphic functions /(^) on M^'^ are as follow 




then in all two cases (i) and (ii), it holds that 



/(^ + e,- + zr,-) = /(^), j = l,---,Ar. 



(3.3) 




2'Ki6jk + 



^k.(^,r) 



% In ■(?(^ + 6^- +zr,-,T) 



(3.4) 
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Differentiating (3.4) witfi respective to again immediately proves tfie formula (3.3) 
for the case (i). The formula (3.4) can be proved for the case (ii) in a similar manner. 
□ 

Theorem 1. Suppose that ■»?(£, e', Ojr) and ■»?(£, e, 0|t) are two Riemann theta 

functions on M^'^, in which e = (si, . . . , En), s' = {e'l, . . . , e'^), and $, = (^i, • • • , ^n), 

= ajX + Ujt + Pj6a + Sj, j = 1,2, ■ ■ ■ , N. Then Hirota bilinear operators D^, -Dj 

and super-Hirota bilinear operator Sx exhibit the following perfect properties when 

they act on a pair of theta functions 

Dx'&{te',0\T)-'&{te,0\T) 

,t^a;Vl^z,q,& - t, < ;||=o v(2^, s' + £, /x/2|2r), 

L M 



J]a,i?(2^,£'-£,-M/2|2T)|4=o 
J]S),^(2£,e'-£,-/u/2|2T)|^=o 



^(2£,e' + e,/u/2|2T), 



(3.6) 



where /u = (//i, • • • ,//Ar), and the notation represents 2^ different transforma- 
tions corresponding to all possible combinations //i = 0, 1; • • • ; //jv = 0, 1. 

In general, for a polynomial operator F{Sx, D^, D^) with respect to S^, and 
Df, we have the following useful formula 



F{Sx,Dx,Dt)i){te\0\T)-i}{te,0\T) 



d{2i,e' + e,n/2\2T), 
(3.7) 

in which, explicitly 

C(£, £', n)= ^ F{M) exp [-27r(r(n - ix/2), n - ix/2) - 2m{n - ix/2, e' - e)\ . 

(3.8) 

where we denote Al = (47ri(n — /u/2, a), AT:i{n — ijl/2,!j:), 47ri(n — /li/2, <t + ^a)). 

Proof. For simplicity we prove the formula (3.6) for one-dimensional case. The 
proof for A^-dimensional case can be performed simply by replacing one-dimensional 
vectors by A'"-dimensional ones. 
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Making use of Proposition 1, we obtain the relation 

A = S,,9{i,e',0\T)-i&{i,s,0\T) 

— ^ 'Dx exp{27rim'(^ + e') — Trm'^r} • exp{27rim(^ + ~ Trm^r}, 

= ^ 27ri(^CT + ^Q;)(m' — m) exp {27ri(m' + m)^ — 27ri(mV + ms) — Trrfm'^ + m^]] 

By shifting sum index as m — I' — m', then 

A= 2%i{a + 9a){2m -l')exp{2%il'^-2%i[m'e' + -m')e]-7rT[m^ + - 

l',m'eZ 

' ^ ^ 27ri(/3a + ^a)(2m' - 2/ - /x) exp{47rie(/ + /^/2) 

/:i=0,l i,m'eZ 

-27ri[mV - (m' - 2Z - //)£] - 7r[m'^ + (m' - 2Z - /x)^} 
Finally letting m' = n + I, we conclude that 



^ 47ri(/3cr + ea)[n - fj,/2] exp{-27rz(n - fJ-/2){e' - e) - 27rr(n - 



^ exp{27ri(Z + /i/2)(2^ + e' + e) - 27rr(Z + /x/2)^ 



^ 2)^i?(2C,e'-e,-///2|2r)|^=o 

M=0,1 



^9(2C,e' + £,M/2|2r), 



by using the following relations 

n + l^{n- pl/2) + {I + ///2), n - I - n ^ {n - ///2) - (Z + n/2). 



In a similar way, we can prove the formula (3.5). The formula (3.7) follows from 
(3.5) and (3.6). □ 

Remcirk 2. The formulae (3.7) and (3.8) show that if the following equations 
are satisfied 

C{s,e',n) = 0, (3.9) 

for all possible combinations /ii — 0, 1; //2 = 0, 1; • • • ; /in — 0, 1, in other word, all 
such combinations are solutions of equation (3.9), then ■i?(^, e', 0|r) and ■j?(^, £, 0|r) 
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are A^-periodic wave solutions of the bilinear equation 

F{St, D„ A)^(£, 0|t) • £, 0|t) = 0. 

We call the formula (3.9) constraint equations, whose number is 2^. This formula 
actually provides us an unified approach to construct multi-periodic wave solutions 
for supersymmetric equations. Once a supersymmetric equation is written bilinear 
forms, then its multi-periodic wave solutions can be directly obtained by solving 
system (3.9). 

Theorem 2. Let C(e,e', /it) and F{Sx, D^, Dt) be given in Theorem 1, and 
make a choice such that e'j — Sj = ±1/2, j = 1, • • • , N. Then 

(i) If F(Sx, Dx, Dt) is an even function in the form 

F{-Sx, -Dx, -Dt) = F{Sx, Dx, Dt), 

then C(e, e', fi) vanishes automatically for the case when X^^^ is an odd number, 
namely 

N 

C(e,e',/u)|^ = 0, for ^/x^ = 1, mod 2. (3.10) 

(ii) If F{Sx, Dx, Dt) is an odd function in the form 

F{-Sx, -Dx, -Dt) = -F{Sx, Dx, Dt), 

then C(e, e', fi) vanishes automatically for the case when is an even number, 
namely 

N 

C{£,£',^i)\^ = 0, for ^/Xj- = 0, mod 2. (3.11) 

Proof. We are going to consider the case where F{Sx, Dx, Dt) is an even function 
and prove the formula (3.9). The formula (3.11) is analogous. Making transforma- 
tion n — — n + fj, {fi = (ni, • • • ,riAr), nj G Z, j — 1, • ■ ■ , N. ), and noting 
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F{Sx, Dx, Dt) is even, we then deduce that 

C(e, e', /x) = F{-M) exp [-27r(T(n - /li/2), n - /u/2) + 27ri(n - /Lt/2, e' - e)] 
— C(e, e', /Lt) exp (47ri(n — /Lt/2, e' — e)) 



/ \ / \ 

= C(e, e', /Lt) exp j ±27ri j exp j ±7ri /i^ j 



= C(£, e', /Lt) exp(±7rz) = — C(s, e', /z), 
which proves the formula (3.10). □ 

Corollary 1. Let e'j — Ej = ±1/2, j — 1, - ■ ■ , N. Assume F{Sx, D^, A) is a 
hnear combination of even and odd functions 

F(5„ Dx, A) = F^{Sx, Dx, A) + F2{Sx, A, A), 

where Fi{Sx, D^, Dt) is even and F2{Sx, D^, Dt) is odd. In addition, C{e,e',ix) 
corresponding (3.8) is given by 

C(e, e', /Lt) = Ci(e, e', /Lt) + C2(e, e', /Lt), 

where 

Ci(e, e', /Lt) = Fi(M) exp [-27r(T(n - /Lt/2), n - /Lt/2) - 27ri(n - /Lt/2, e' - e)] , 



C2(e, e', /Lt) = F2(M) exp [-27r(T(n - /Lt/2), n - /Lt/2) - 27ri(n - /Lt/2, e' - e)] . 
Then 



AT 



C(e,e',/Lt) = C2(e,e',/Lt) for ^/x^- = 1, mod 2, (3.12) 

N 

C(e,e',/Lt) = Ci(e,e',/Lt), for ^lij = 0, mod 2. (3.13) 

Proof. In a similar to the proof of Theorem 2, shifting sum index as n = — n + /Lt, 

and using Fi(5'x, A, A) even and F2{Sx, A, A) odd, we have 
C(e, e', /Lt) = Ci(e, e', /Lt) + C2(e, e', /Lt) 

= [Ci(£,e',/u) - C2(e,e',/ti)]exp ^±7ri^/Xjj . ^^'^^"^ 
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Then for X^jLi/^j — 1; mod 2, the equation (3.15) gives 

which imphes the formula (3.12). The formula (3.13) is analogous. □ 

The theorem 2 and corollary 1 are very useful to deal with coupled super-Hirota's 
bilinear equations, which will be seen in the following section 5. 

By introducing differential operators 

V = (%,%,■■■ ,9^^), 

(9^ = ai% + 025^2 H h onO^^ = ck ■ V, 

dt = + u;2% H \- (^Nd^N = ^ • V, 

2) = ((Ti + ^ai)% + (as + ^02)% + ■ ■ ■ + ((T^ + ^«7v)5civ = (o" + ^ck) ■ V, 
then we have 

T>d',dmt r) = [{a- + 0a) ■ V](« ■ V)^cu ■ VYHt r) 

= (cr ■ V)(a ■ V)'^(lc; • V)%(^, r) + ^(a ■ Vf^^uj ■ V)%(^, r), (3.15) 
j,/c = 0,l,--- . 

4. = 1 supersymmetric Sawada-Kotera-Ramani equation 

The supersymmetric Sawada-Kotera-Ramani equation takes the form 

(j)t + [1O(D0)S)V + 5(DV)0 + 15(D0)V] + = 0, (4.1) 

where = 0(x, t, 6) : M^'^ — t- M^'^ is fermionic superfield depending on usual inde- 
pendent variable x, t and Grassmann variable 6. This equation was first proposed 
by Carstea [12] • The soliton solutions, Lax representation and infinite conserved 
quantities of the equation have been further obtained recently [55t [56] . Here we are 
interested in quasi-periodic wave solutions to the supersymmetric equation (4.1). 
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We will show that the soliton solutions can be obtained as limiting case of these 
quasi-periodic solutions. 

To apply the Hirota bilinear method in superspace for constructing multi-periodic 
wave solutions of the equation (4.1), we hope more fee variables and consider a gen- 
eral variable transformation 

<l>^<l>o + 2^Hnf{x,t,e), (4.2) 

where f{x,t,6) : M^'^ — )> M.]f is an even superfield and 0o = 0o(^) • H^a"*^ ^^a^ 
is an odd special solution of the equation (4.1). Substituting (4.2) into (4.1) and 
integrating with respect to x, we then get the following super Hirota's bilinear form 

F{S,, D,, Dt)f ■ f = {S,Dt + S^Dl + IbcgDl + bct>oDl + c)f-f = 0, (4.3) 

where c = c{9, t) : M.^^ — >■ ]R°'^ is an odd integration constant. In the special case 
when 00 = c = 0, starting from the bilinear equation (4.3), it is easy to find that the 
equation (4.1) admits one-soliton solution (also called one-supersoliton solution) in 
superspace M^'^ over two-dimensional Grassmann algebra Gi{a) 

(j)i = 2^Hn{l + e"), (4.4) 

with phase variable rj — px — pH -\- q9a + r with p,q,r e Aq. While two-soliton 
solution (super two-soliton solution) reads 

02 = 2S)^ln(l + e*?! + e''^ + e'^'+'^'+^''), (4.5) 

with rjj = pjX — p^t + QjOa + rj, j = 1,2 and 

gAi2 ^ jPi - P2)^{pI - P1P2 + pI) a ^ ^^^ Piq2 - P2qi \ 
{pi+P2)'^{pi+PiP2 + pl) \ P1-P2 /' 

and here pj, qj, rj & Aq, j — 1,2 are free constants. 

Next, we turn to see the periodicity of the solution (4.2), the function / is chosen 
to be a Riemann theta function, namely, 

f{x,t,9)=i^{tr). 
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where phase variable ^ is taken as the form ^ — {Ci, ■ ' ' ■> Cn)'^, Cj — + ^j't + 
/3j9a + 5j, j — 1,2, - ■ ■ ,N. With Proposition 3, we refer to 

TV 

= 00 + 2 5] akiPia + 9ai)dl^^ r), 

k,l=l 

which shows that the solution is indeed a quasi-periodic function with 2N funda- 
mental periods {e^-, j = 1, • • • , A^} and {iTj, j — 1, - ■ ■ , N}. The quasi-periodic 
means that (f) is periodic in each of the N phases {^j, j — 1, - ■ ■ , N}, if the other 
N — 1 phases are fixed. 

3.1. One-periodic waves and asymptotic analysis 

We first consider one-periodic wave solutions of the equation (4.1). As a simple 
case of the theta function (3.1) when = 1, s = £ = 0, we take / as 

f{x,t,0)^ ^iC, r) = Y, exp(27rme - TrnV) , (4.6) 

neZ 

where the phase variable ^ — ax + out + f39a + S, and the parameter r > 0. 

Next, we let the Riemann theta function (4.6) be a solution of the bihnear equa- 
tion (4.3). By using Theorem 1 and the formula (3.9), the following two equations 
( corresponding to = and 1 respectively) should be satisfied 

^{-(47r(n - ii/2)f{^(7 + ea)uj - (47r(n - n/2)f{^(T + ea)a^ - 15(47r(n - ii/2)fa^(t)l 

n&L 

+ 5(47r(n - /x/2)) V^o + c} exp(-27rT(n - Ai/2)^) = 0, = 0, 1. 

(4.7) 

We introduce the notations by 

A = e--l\ ^i(e. A) = ^(2e, 0, 0|2t) = A^"' exp(4z7rnO, 
^2(e, A) = 0, -1/2|2t) = A(2"-^)' exp[2z7r(2n - 1)^], 
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the equation (4.7) can be written as a linear system about cu, c 

{/3a + ea)d'-uj + (/3(7 + ea)a^'df + ^a^'d'-^l + Sa^'i?^ Vo + ^jC = 0, j = 1, 2, 

(4.8) 

where a; e Aq is even and c, 00 G Ai are odd, and we have denoted the derivative 
value of 'dj{^i A) at ^ = by simple notations 

^,; = ^,;(o,A) = !^|,.„. j = i,2. 

Moreover, we see that the functions and their derivatives are independent of 
Grassmann variable 6 and anticommuting number a. 

We take 0o = for the simplicity. It is obvious that the coefficient determinant 
of the system (4.8) is nonzero and {aHf\aHfY ^ 0, therefore the system (4.8) 
admits a solution 

where cu is independent of Grassmann variable 9 and auticommuting number a, and 
parameter a is free. 

In this way, a one-periodic wave solution of the equation (4.1) is explicitly ob- 
tained by 

= 2D-Mni?({,r), (4.10) 

with the theta function ■j?(^,t) given by (4.6) and parameters u, c by (4.9), while 
other parameters a, /3,t,S G Aq are free. Among them, the three parameters a and 
r completely dominate a one-periodic wave. 

In summary, one-periodic wave (4.10) possesses the following features: 

(i) It is one-dimensional, i.e. there is a single phase variable ^. Moreover, it has 
two fundamental periods 1 and ir in phase variable ^, but it need not to be periodic 
in X, t and 9 directions. 
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Figure 1. A one-periodic wave for M = \ supersymmetric Sawada-Kotera- 
Ramani equation with parameters: a — 0.1, r — 32, cr = 0.013. (a) Perspective 
view of wave, (b) Overhead view of wave, with contour plot shown. 

(ii) It can be viewed as a parallel superposition of overlapping one-soliton waves, 
placed one period apart ( see (a) and (6) in Figure 1 ). 

(iii) Different form the purely bosonic case, it is observed shows that there is an 
influencing band among the one-periodic waves under the presence of the Grassmann 
variable (in contour plot, the bright hexagons are crests and the dark hexagons are 
troughs). The one-periodic waves are symmetric about the band but collapse along 
with the band. Furthermore, the amplitudes of the quasi-periodic waves increase as 
the waves move away from the band ( see (a) and (6) in Figure 1 ). 

(iv) The quasi-periodic wave will degenerate to pure bosonic quasi-periodic wave 
when Q becomes small ( see Figure 2 ). 

In the following, we further consider asymptotic properties of the one-periodic 
wave solution. Interestingly, the relation between the one-periodic wave solution 
(4.10) and the one-soliton solution (4.4) can be established as follows. 

Theorem 3. Suppose that w and c are given by (4.9), and for the one-periodic 
wave solution (4.10), we let 




r -|- TTT 



(4.11) 



a = 



2Txi 
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Figure 2. A purely bosonic case of one-periodic wave to A/" = 1 supersym- 
metric Sawada-Kotera-Ramani equation with parameters: a = 0.1, r — 32, a — 
0.013, 9 ~ 0.00000001. (a) Perspective view of wave, (b) Overhead view of wave, 
with contour plot shown. 

where the p,q and r are given in (4.4). Then we have the following asymptotic 
properties 

c — > 0, 2niC, — nr — )■ t] = px — p^t + qOa + r, 

??(^,r) — ^ 1 + 6^ as A ^ 0. 
In other words, the one-periodic solution (4.10) tends to the soliton solution (4.4) 

under a small amplitude limit, that is, 

— ^01, as A -> 0. (4.12) 

Proof. Here we will directly use the system (4.8) to analyze asymptotic prop- 
erties of one-periodic solution (4.10), which is more simple and effective than our 
original method by solving the system Since the coefficients of system 

(4.8) are power series about A, its solution (w, c) also should be a series about A. 
We explicitly expand the coefficients of system (4.8) as follows 
t9i = 1 + 2A^ + ■ ■ ■ , d'[ = -327r2A^ + ■ ■ ■ , 

^^S^) = -81927r^A^ + ■ ■ ■ , d2 = 2X + 2\^ + ■■■ (4.13) 

^9^' = -87r2A + • ■ ■ , df^ = -12Stx^\ + --- . 
Let the solution of the system (4.8) be of the form 

uj = bjQ + uJiX + UJ2\^ H = Wo + o(A), 

(4.14) 

c = Co + ciA + C2A H = Co + o(A). 
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Substituting the expansions (4.13) and (4.14) into the system (4.8) (the second 
equation is divided by A ) and letting A — > 0, we immediately obtain the following 
relations 

Co = 0, -87r^(^(7 + ea)uJo + 2co - 1287r^(^(7 + ea)a^ = 0, 
which has a solution 

Co = 0, wo = -IGTr^a^ (4.15) 
Combining (4.14) and (4.15) then yields 

c — y 0, 2Tiiuj — y -32in^a^ = -p^, as A ^- 0. 

Hence we conclude 

^ = 2m^ — TTT — px + 2'!riujt + qOa + r 

(4.16) 

— > px — p t + qOa + r = 7], as A ^ 0. 

It remains to consider asymptotic properties of the one-periodic wave solution 
(4.10) under the hmit A — > 0. By expanding the Riemann theta function ■j?(^, r) and 

by using (4.16), it follows that 

^?(^, r) = 1 + X\e^^'^ + e-2"'«) + A^(e^"'« + e-^"'«) + • • • 

= 1 + e^" + A^(e-« + e^^) + X'\e-^^ + e'^) + ■■■ 

— ^ 1 + e^' — yl + e", as A ^ 0, 
which together with (4.10) lead to (4.12). Therefore we conclude that the one- 
periodic solution (4.10) just goes to the one-soliton solution (4.4) as the amplitude 
A^ 0. □ 

3.2. Two-periodic wave solutions and asymptotic analysis 

We proceed to the construction of the two-periodic wave solutions to the su- 
persymmetric Sawada-Kotera-Ramani equation (4.1), which are a two-dimensional 
generalization of one-periodic wave solutions. The two-periodic waves of interest 
here have three-dimensional velocity fields and two-dimensional surface patterns. 
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For the case when A^ = 2,s = e = 0in the Riemann theta function (3.1), we 
takes / as 

/ = ^) = E exp{27rz(^, n) - 7r(rn, n)}, (4.17) 

where n = (ni, ^2)^ G Z^, £ = (^1, ^2)^ e = Q;jX + a;ji + ^j6l(7 + 5j-, j = 1, 2; 

The matrix r is a positive definite and real- valued symmetric 2x2 matrix, which 
can takes of the form 

T = (Tij)2x2, T12 = T21, Til > 0, T22 > 0, T11T22 - ri2 > 0- 

Next, we explore the conditions to make the Riemann theta function (4.17) 
satisfy the bilinear equation (4.3). Theorem 1 and the formula (3.9) give rise to the 
following four constraint equations 

^ [-An^{2n - n,cr + 9a){2n - /u, u>) - 647r^(2n - /x, af{2n - tJ,,cr + 9ct) 

ni,n2& 

- 607r^(2n - /j,, a^cpl + 807r^(2n - /j,, a) Vo + c] exp{-27r(r(n - n/2),n - n/2)} = 0, 

(4.18) 

where fi = (/^i, 112) takes all possible combinations of /^i, /i2 = 0, 1. 
By introducing the notations 
Xki = e-'^"'='/^ k,l ^ 1,2, (All, A12, A22) 

A) = m, 0, -3,/2|2r) = Yl exp{47rz(t n - Sj/2)} fl Ag"'=-^^'=)(^"^-^-.'), 

ni,n2&Z k,l=l 

Sj^(sj,i,Sj,2), 1,2,3, A, si = (0,0), S2 = (l,0), S3 = (0,1), S4 = (l,l), 

then by using (3.15), the system (4.18) can be written as a linear system 
[(/3(7 + ^a) • V](w • V)^?, (0, A) + [(/3a + ^a) • V](a • V)%(0, A) 

(4.19) 

+ 15(a • V)%(0, A)0g + 5(a • V)%(0, A)(/)o + ^j(0, A)c = 0. 
This system is easy to be solved in such a way: by solving a quadratic equation 
with one unknown; a;i,a;2 and c by solving a linear system. With such a solution 
(a;i,a;2, 0o, c), we then get an exact two-periodic wave solution 

(/) = 0o + 2S)'ln^?(tr), (4.20) 
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Figure 3. A two-super periodic wave for TV = 1 supersymmetric Sawada- 
Kotera-Ramani equation with parameters: ai — 0.1, a2 — —0.1, th — 2,ti2 = 
0.2, T22 ~2,ai = 0.1, CT2 = 0.1. (a) Perspective view of wave, (b) Overhead view 
of wave, with contour plot shown. 



with t) and ijji,co2,(po,c given by (4.17) and (4.19), respectively, while other 
parameters /32, Oi, a2, Tii, T22, T12, ^1,^2 e Aq are free. 

In summary, two-periodic wave (4.20), which is a direct generalization of one- 
periodic wave, has the following features: 

(i) The two-periodic wave solution is genuinely two-dimensional. Its surface 
pattern is two-dimensional, namely, there are two phase variables ^1 and ^2- 

(ii) It has two independent spatial periods in two independent horizontal direc- 
tions. It has 4 fundamental periods {61,62} and {iTi,iT2} in (^1,^2)- It is spatially 
periodic in two directions ^1,^2, but it does not need periodic in the all x-, t- and 
^-directions. 

(iii) As in the case of on-periodic waves, there is an influencing band among the 
two-periodic waves under the presence of the Grassmann variable. ( see Figure 3 ). 

At last, we consider the asymptotic properties of the two-periodic solution (4.20). 
In a similar way to Theorem 3, we can establish the relation between the two-periodic 
solution (4.20) and the two-soliton solution as follows. 
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Theorem 4. Assume that {uji,uj2, 0o, c) is a solution of the system (4.19). We 
choose the parameters in the two-periodic wave solution (4.20) as follows 



with the Pj,qj,rj,j = 1, 2 as those given in (4.5). Then under constraint ai(32 — 
q;2/3i, we have the following asymptotic relations 

00 — > 0, c — > 0, — TTTjj — > pjX — p^t + QjOa + Tj ^rjj, j = 1, 2, 

^{^, r) ^ 1 + + 6"^ + e''i+^2+Ai2^ ^ ^^^^ _^ q 

(4.22) 

So the two-periodic wave solution (5.20) tends to the two-soliton solution (4.5), 
namely, 

— > 02, as All, A22 0. 

Proof. Prom (4.21), the constraint q;i^2 = Q;2/5i leads to piq2 — P2Q1 — 0, which 
implies that T12 = —^12/271 is independent of Grassmann variable 9 according to 
(4.5). 

In the same way as the proof of Theorem 3, we expand the Riemann function 
^(Ci> C2, t) in the following form 

i?(£, r) = 1 + (e^'^'^i + e-2'^^^i)e-''^" + {e^'''^^ + g-^'^^^^^g-TrT^a 

_^^g27ri(a+6) _^ g-27ri(6+6))g-T(Tll+2Tl2+T22) _^ 

Further by using (4.21) and making a transformation a)j = 2mujj,j — 1, 2, we get 
^?(£, r) = 1 + e«"i + + eSi+^^-^-^i^ + ^4^g-li ^ ^4^g-6 ^ A^iA^ae"^^"^''"'^^'' + 
^ 1 + e^'i + e^"^ + e^'^+^"^+^i% as Aii,A22^0, 
where — pjX + oujt + Pj9a + rj, j — 1,2. 
It remains to prove that 

c — >0, Uj -p], ij — > rij, j = 1, 2, as An, A22 ^ 0. (4.23) 
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As in the case when = 1, we let the solution of the system (4.19) be the form 

^1 = ^1,0 + 1^1,1-^11 + '^1,2-^22 + ■ ■ ■ = Wi^o + o(Aii, A22), 
UJ2 = ^2,0 + 1^2,1-^11 + lx'2.2A22 + ■ ■ ■ = U!2,0 + o(Aii, A22)5 

(4.24) 

00 = 00,0 + 0O,lAll + 0O,2A22 H = 00,0 + o(Aii, A22), 

c = Co + ciAii + C2A22 H = Co + o(Aii, A22)- 

Expanding functions -(9^, j = 1, 2, 3, 4 in the system (4.19), together with substitution 
of assumption (4.24), the second and third equation is divided by An and A22, 
respectively; the fourth equation is divided by A11A22, and letting Aii,A22 — > , 
we then obtain 

Co = 0, -87r(/3i(T + eai)uji - 1287r^a^(/3i(T + Ooi) - 12O7ra^0^o + 20o,o = 0, 

- 87r(/32fT + ^^2)^2 - 1287r^a^(/32a + 602) - UOiralcPlo + 20o,o = 0, 
which has solution 

Co = 0, 00,0 = 0, ujifi = -Wn'^al, a;2,o = -Wn^a^. (4.25) 

The expressions (4.24) and (4.25) implies that 

00 = o(Aii, A22) — >0, c = o(Aii,A22) — >0, = -levr^ai + o(Aii, A22) — > -Wn 
W2 = -16vr'^a2 + ^^(-^11) -^22) — > -IQir'^al, as Aii,A22-^0, 

thus proving (4.23). We conclude that the two-periodic wave solution (4.20) tends 
to the two-soliton solution (4.5) as An, A22 — )■ 0. □ 

4. J\f = 2 supersymmetric KdV equation 

We consider M = 2 supersymmetric KdV equation 

<Pt = -0XXX + 3(0DiD20). + ^(a - l)(S)iD20').. + 3a0202, (5.I) 

which was originally introduced by Laberge and Mathieu [571 EH] • In the equation 
(5.1), = (f){x, t, 61, 62) : M^'^ — )• R^''^ is a superboson function depending on tempo- 
ral variable t, spatial variable x and its fermionic counterparts 61, 62. The operators 
Di and 'D2 are the super derivatives defined by S)i = de-^ + Oid^, 2)2 = dg^ + 620^ 
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and a is a parameter. The equation (5.1) is called supersymmetric KdVa equation 
[5U] . For the cases when a = 1 and a = 4, the Lax representation, Hamiltonian 
structure, Painleve analysis and soliton solutions of the equation (5.1) can refer to, 
for instance, papers [S7]-[5S]. 

Here we are interested in quasi-periodic wave solutions to the supersymmetric 
equation (5.1) by using Theorem 1 and 5. We only consider the case when a = 1, 
so the equation (5.1) reduces to 

0t = + 3(0S)i2)20). + 30202_ (5 2) 

To apply the Hirota bilinear method for constructing multi-periodic wave so- 
lutions of the equation (5.2), we add two variables and consider a general variable 
transformation 

/ 

(^ = u + 92V, u = idxln—, v = vo — Ox's lTi(f g) , (5.3) 

9 

where u{x,t,6i), f{x,t,6i), g{x,t,6i) : M^'"*^ — )■ and v{x,t,9i),Vo = Vo{6i) : 

M^'^ —7- R°'^ is a special solution of the equation (5.2). Hereafter we use D = Di for 
simplicity. Substituting (5.3) into (5.2), we then get the following bilinear form 
FiDt,D.x)f-g = iDt + Dl)f-g = 0, 

(5.4) 

GiSx, Dt, a,)f ■ g = {S,Dt + S^Dl + 3voDl + c)f-g = 0, 

where c = c{6i,t) : R^'"*^ — > ]R°'^ is an odd integration constant to variable x; The 
equation (5.4) is a type of coupled bilinear equations, which is more difficult to be 
dealt with than the bilinear equation (4.3) due to appearance of two functions and 
two equations. We will take full advantages of Theorem 2 to reduce the number of 
constraint equations. 

Now we take into account the periodicity of the solution (5.3), in which we take 
/ and g as 

/ = ^(^ + e,r), g = i^i^ + h,T), e,/iGZ^, 
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where phase variable ^ is taken as the form ^ = (^i, ■ ■ ■ ,^n)'^, = + ^jt + 
(3j9ia + 6j, j = 1,2, ■ ■ ■ , N. By means of Proposition 3, we deduce that 

N 

u + 62V, u = i 



> a/cC'fi. In — — ; r 

^ m + h,T) 



N 

v = v,~2Y, + dai)dl^^ + e, t)^(^ + r)], 

which indicates that the solution is a quasi-periodic function with 2N fundamental 
periods {sj, j = 1, - ■ ■ , N} and {iTj, j = 1, ■ ■ ■ , A^}. 

In the special case of = c = 0, starting from the bilinear equation (5.4), Zhang 
et al. found that the equation (5.2) admits one-soliton solution [59] 

01 = id^ In ^ + e^lv^ - Hh9i)l (5.5) 
9i 

with 

/l = l + e^ gi = l-e^ 

and phase variable r] = px — pH + qOia + r with p,q,r G Aq. While two-soliton 
solution takes the form 

01 = id. In ^ + ^2 bo - 5,2) ln(/2^2)], (5.6) 

^2 



with 

/2 = 1 + + 6"^ + e''^+''^+^i% 
and ?7j = PjX — pp + qjOia + r^, j = 1, 2, 



(5.7) 



(pi+P2)^ {pi+P2y 

here pj, g^, G Aq, j = 1, 2 are free constants. 

4.1. One-periodic waves and asymptotic analysis 
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We first construct one-periodic wave solutions of the equation (5.2). As a simple 
case of the theta function (3.2) when N — l,s — 0, we take / and g as 

/ = ^?(e, 0, 0|t) = exp(27rme - TrnV), 

g = ^(e, 1/2, 0|t) = ^ exp(27rzn(e + 1/2) - TrnV) (5 3) 

= y^(-l)"exp(27rw^ - Trn^r), 

nGZ 

where the phase variable ^ — ax + out + pOia + S, and the parameter r > 0. 

Due to the fact that F{Dt, D^) is an odd function, its constraint equations in the 
formula (3.10) vanish automatically for = 0. Similarly the constraint equations 
associated with G{Sx, Dt, Dj.) also vanish automatically for /j, — 1. Therefore, the 
Riemann theta function (5.8) is a solution of the bilinear equation (5.4), provided 
the following equations 

^{27ri(2n - //)a; - i{27raf{2n - nf} exp[-27rT(n - ///2)^ + 7ri{n - ///2)] |^=i = 0, 
^{-[27r(2n - fJ,)]'^{pa + 0ia)uj + {2n{2n - fJ,)]'^{Pa + eia)a^ - {2n{2n - fJ,)a\'^vo + c} 
X exp(-27rT(n - i^/2)^ + 7ri{n - /x/2))|^=o = 0. 

(5.9) 

We introduce the notations by 

MC, A) = m, 1/4, -1/2|2t) = J2 A^'"~'^' exp[4m(n - l/2)(e + 1/4)], 
M^, A) = m, 1/4, 0|2t) = Yl exp[4mn(e + 1/4)], 

the equation (5.9) can be written as a linear system about ou, c 
^[u + a^'l' = 0, 

(5-10) 

{/3a + 9ia)&2i^ + + {/3a + 9ia)a^^r + a'^^2Vo = 0. 



^'. = ^'.(0,X) = ^-^U^, J = 1,2. 
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where a; e Aq is even and c, vq G Ai are odd, and we have denoted the derivative 
value of A) at ^ = by simple notations 

Moreover, we see that the functions 1?^ and their derivatives are independent of 
Grassmann variable 6 and anticommuting number a. 

We take vq — ^a{a + 9a), 7 e Aq for the simplicity, then the system (5.10) 
admits a solution 

where cu is independent of Grassmann variable 9 and anticommuting number a. In 
this way, a one-periodic wave solution reads 

= id. In ^^^'^ifl^^ + ^2{^o - 9,2) ln[^(^, 0, 0|t)^(^, 1/2, 0|r)]}, (5.12) 

where parameters uj and c are given by (5.11), while other parameters a, (3,T,d e Aq 
are free. Among them, the three parameters a and r completely dominate a one- 
periodic wave. 

In summary, one-periodic wave (5.12) has the following features: 

(i) It is onc-dimciisional and has two fundamental periods 1 and ir in phase 
variable ^. It can be viewed as a parallel superposition of overlapping one-soliton 
waves, placed one period apart (see Figure 5-7). 

(ii) As in the case of the supersymmetric Sawada-Kotera-Ramani equation, there 
is also an influencing band among the real part of one-periodic waves for the super- 
symmetric KdV equation under the presence of the Grassmann variable (see Figure 

4). 

(iii) It was not observed that influencing band appears among the imaginary 
part and modulus of the one-periodic wave. Moreover, they seem to have the same 
shape from the observation of their plots (see Figures 5 and 6) . 
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(a) (6) 




Figure 4. Real part of one-periodic wave for A/" = 2 supersymmetric KdV 
equation with parameters: a = 0.1, t = 3, ui — 0.01. (a) Perspective view of 
wave, (b) Overhead view of wave, with contour plot shown. 

(a) (6) 




10 1 



Figure 5. imaginary part of one-periodic wave for J\f = 2 supersymmetric 
KdV equation with parameters: a ~ 0.1, t — 3,ai = 0.01. (a) Perspective view 
of wave, (b) Overhead view of wave, with contour plot shown. 

(a) (6) 




Figure 6. Modulus of one-periodic wave for J\f = 2 supersymmetric KdV 
equation with parameters: a — 0.1, t — 3,ai = 0.01. (a) Perspective view of 
wave, (b) Overhead view of wave, with contour plot shown. 
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In the following, we further consider asymptotic properties of the one-periodic 
wave solution. The relation between the one-periodic wave solution (5.12) and the 
one-sohton solution (5.5) can be estabhshed as follows. 

Theorem 5. Suppose that u and c are given by (5.11). In the one-periodic 
wave solution (5.12), we choose parameters as 

7 = 0, «=— , /3 = — . 5= , (5.13 

2m Zni 2m 

where the p, g and r are the same as those in (5.5). Then we have the following 
asymptotic properties 

c^O, i^1±^, /^l + e^ g^l-e\ as A^O. 

In other words, the one-periodic solution (5.12) tends to the one-soliton solution 
(5.5) under a small amplitude limit , that is, 

— as A^^O. (5.14) 

Proof. Here we will directly use the system (5.10) to analyze asymptotic proper- 
ties of one-periodic solution (5.12). We explicitly expand the coefficients of system 
(5.10) as follows 

= -47rA + 127rA^ + • • • , = IGtt^A + 4327r=^A'^ + • • • , 

= -2 + 2A^ + --- , i?^' = 327r'A^ + • • • , (5.15) 
^""^ = 5127r^A^ + • • • 
Suppose that the solution of the system (5.10) is of the form 

u ^uo^- ui\ U2\^ H = u;o o(A), 

(5.16) 

c= cq + ciA -\- C2A H = Co o(A). 

Substituting the expansions (5.15) and (5.16) into the system (5.11) and letting 
A — > 0, we immediately obtain the following relations 

- A-KUjQ + 167r^Q;^ = 0, cq = 0, 
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which has a solution 

Co = 0, Wo = ^-K^a^. (5.17) 
Combining (5.13) and (5.17) leads to 

c — 0, 2muj — > Sn^ia^ = -p^, as A ^- 0, 

or equivalently 

^ = 2m^ — TTT = px + 27iiuJt + q9ia + r 

(5.18) 

— > px — p t + q9\a + r — T], as A ^ 0. 

It remains to identify that the one-periodic wave (5.12) possesses the same form 
with the one-soliton solution (5.5) under the limit A ^ 0. For this purpose, we start 
to expand the functions / and g in the form 

/ = 1 + A2(e2"« + e-2"'«) + A^(e^'^'^ + e'^^'^) + ■■■ . 

g^l- A'(e''^'^ + e-^""'^) + X^e^""'^ + e-^"'«) + • • • . 

By using (5.13) and (5.17), it follows that 

f = l + J + X\e-^ + e^^') + X^^e-^^ + e^^) + ■■■ 

— >l + e^ — >l + e'^, as A ^ 0; 

, , (5.19) 
g^l-e^ + X\e'^ - e-«) + X'^ie''^ - e'^) + ■■■ 

— ^ 1 - — ^ 1 - as A ^ 0. 
The expression (5.13) follows from (5.19), and thus we conclude that the one-periodic 
solution (5.12) just goes to the one-soliton solution (5.5) as the amplitude A — > 0. 
□ 

4.2. Two-periodic waves and asymptotic properties 



We now consider two-periodic wave solutions to the supersymmetric KdV equa- 
tion (5.2). For the case when = 2, s = 0, e = 1/2 = (1/2, 1/2) in the Riemann 
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theta function (3.2), we choose / and g to be 

/ = ^{t 0, 0|t) - exp{27rz(^, n) - 7r(Tn, n)}, 

^ = d{t 1/2, 0|t) = exp{27ri(^ + 1/2, n) - 7r(Tn, n)} (5 20) 

= Y (-1)'''+''' exp{27ri(^, n) - n{Tn, n)} 

where we denote n — (711,^2) G Z"^, ^ = (^1, ^2) € C^, = ajX + 0;^^ + /5j^icr + 
5j, j — 1,2, and a = (q;i,q;2), u? = {uji,uj2), (3 — {^1,^2) £ C^,; The matrix r is a 
positive definite and real- valued symmetric 2x2 matrix, which can take the form 

T = (Tii)2x2, T12 = r2i. Til > 0, r22 > 0, T11T22 " rf^ > 0. 



According to Theorem 5, constraint equations associated with F{Dt, D^) vanish 
automatically for (111,112) — (0,0), (1,1), and the constraint equations associated 
with G{Sx, Df, D^) vanish automatically for (/ii,/U2) = (1,0), (0, 1). Hence, making 
the theta functions / and g satisfy the bilinear equation (5.4) gives to the following 
constraint equations 

Y [27ri(2n - /u, u?) - 87r^i(2n - /u, af] exp{-27r(T(n - /u/2), n - /u/2) 

2 

+ TTi ^(rij - /ij/2)}|^=(^i,^2) = 0, for (//i, ^2) = (0, 1), (1, 0). 

(5.21) 

and 

Y [-47r^(2n - n,a/3 + 9iOt){2n - fi,u>) + 167r^(2n -fi,af3 + di(x){2n - n, 

- 47r^(2n - /jl, a^vo + c] exp{-27r(T(n - n/2),n - /Lt/2) 

2 

+ ni Yi'^j - /^j/2)}U=(w,/^2) = 0, for (/xi, /X2) = (0, 0), (1, 1). 

(5.22) 
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Next, introducing the following notations 



Xki = e-'^"'='/^ k,l ^ 1,2, \^ (An, A12, A22) 

2 

^?,(tA)=^(2e,l/4,-s,/2|2r)= J] exp{Am{^ + l/4,n - sj/2)} JJ X^^^^-'^ 

n\,n2&Z k,l=l 

Si = Ki,s,-,2), i = 1,2,3,4, si = (0,l), S2 = (l,0), S3 = (0,0), 54 = (1,1), 



then by using (3.15), the system (5.21) and (5.22) can be written as a linear system 
(a; ■ V)^9,(0, A) + (a ■ V)%{0, A) = 0, j = 1, 2, 

[(a/3 + ^la) • V](u; • V)i?,(0, A) + [(a/3 + ^la) • V](a • V)'^?, (0, A) (5.23) 

+ (a • V)'^?,(0, A)t;o + 1?,(0, A)c = 0, j = 3, 4. 

This system can be solved in such a way: After we obtain u,i , U2 form the first two 
equations. We substitute them into last two equations to get Vq, c. With the solution 
{uji,uj2,vq,c), we get a two-periodic wave solution to the supersymmetric equation 
(5.2) 

= afi^; ^o°LT^ ^ + ^2{^o - 9,2) ln[^(£, 0, 0|t)^(£, 1/2, 0|t)]}, (5.24) 

^^(€, 1/^, ufj 

where parameters a;2, vo and c are given by (5.22), while other parameters ai, ui, ol\, 
0L2-, Til, T22, T12, ?>\ and (^2 are free. 

In summary, the two-periodic wave (5.24), which is a direct generalization of 
one-periodic waves, has the following features: 

(i) Its surface pattern is two-dimensional, namely, there arc two phase variables 
^1 and ^2- It has 4 fundamental periods {61,62} and {iri,ir2} in (^1,^2), and is 
spatially periodic in two directions ^1, ^2- Its real part is not periodic in Q direction, 
while its real part, imaginary part and modulus are all periodic in x and t directions. 

(iii) There is also an influencing band among the Real part of two-periodic waves 
for the supersymmetric KdV equation under the presence of the Grassmann variable 
( see Figure 7 ). 
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(a) (6) 




Figure 7. Real part of two-periodic wave for J\f — 2 supersymmetric KdV 
equation with parameters: ai = 0.1, Q!2 = 0.2 th = 3, T12 = 0.2, T22 — S,(7i = 
0.2, f72 = —0.1. (a) Perspective view of wave, (b) Overhead view of wave, with 
contour plot shown. 



(a) (6) 
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Figure 8. imaginary part of two-periodic wave for Af ^ 2 supersymmetric 
KdV equation with parameters: ai — 0.1, a2 — 0.2 th — 3, T12 = 0.2, T22 = 
3,(Ti = 0.2, (72 = —0.1. (a) Perspective view of wave, (b) Overhead view of wave, 
with contour plot shown. 

(iv) It was not found that influencing band appears among the imaginary part 
and modulus of two-periodic waves to the supersymmetric KdV equation ( see Fig- 
ures 8 and 9 ). 

Finally, we consider the asymptotic properties of the two-periodic solution (5.24). 
In a similar way to Theorem 5, we can establish the relation between the two-periodic 
solution (5.24) and the two-soliton solution (5.6) as follows. 



(a) (6) 
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Figure 9. Modulus of two-periodic wave for A/" = 2 supersymmetric KdV 
equation with parameters: ai = 0.1, a2 — 0.2 th — 3, T12 = 0.2, T22 = 3,cti = 
0.2, (72 = —0.1. (a) Perspective view of wave, (b) Overhead view of wave, with 
contour plot shown. 



Theorem 6. Assume that {u)i,U2,Vo,c) is a solution of the system (5.23). In 
the two-periodic wave solution (5.24), a choice of parameters is given by 

= /?, = ^. 6, = ^^, = J = 1,2, (5.25) 

2m 2m 2m 2tc 

with the Pj,qj,rj G Ao,j = 1,2 and A12 as those given in (5.6). Then under the 
constraint ai/32 = «2/3i, we have the following asymptotic relations 

vo^O, c^O, 0^ \ J = 1,2, 

2tii 

f — >l + e"' +6"'+ e''i+'?2+^i2^ g — ^ 1 _ e''^ - 6"^ + e^'i+^^+^i^ (5.26) 
as All, A22 0. 

So the two-periodic wave solution (5.24) just tends to the two-soliton solution (5.6) 
under a certain limit 

(f) — )■ 02, as All, A22 0. 



Proof. From (5.25), the constraint aif32 = 012^1 leads to piq2 — P2(li = 0, which 
implies that ri2 = — Ai2/27r is independent of Grassmann variable 6 according to 
(5.7). 
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Using (5.20), we explicitly expand the functions / and g in the following form 

_^(g27ri(ei+6) _|_ g-27ri(a+6)-)g-7r(Tii+2Ti2+r22) _^ 

_^(g27ri(6+6) _^ g-27ri(6+6))g-^(ni+2ri2+r22) _^ 

Further using (4.5) and making a transformation Uj = 2'n'iuj,j = 1, 2, we infer that 

/ = 1 + + + e«^+«"^-2^^i2 + A'^^e-^^ + A^ae"^' + XtAe'^'-^''^""^'' + ■■■ 

— > 1 + e^"i + e^"^ + e^'i+^'^+^i^ as An, A22 ^ 0, 
^ = 1 - e^'i - e^"^ + e^'i+^'^-'^^^i^ - AfiC-^'i - A^ae"^"' + A^iA^se"^"'"^"'"'''^'' + • • • 
1 - - e«'^ + e€i+«2+^i^ as An, A22 ^ 0, 

where — pjX + oujt + qjOia + rj, j — 1,2. 
It remains to prove that 

c — )■ 0, Uj — y ij — > Tjj, j = 1, 2, as An, A22 ^ 0. (5.27) 

As in the case when A?" = 1, we let the solution of the system (5.23) be the form 

^1 = ^1,0 + ^1,1^11 + U;2,2A22 + o(Aii, A22), 
UJ2 — UJ2,0 + ^2,1 All + ^'-'2,2^22 + o(Aii, A22), 

(5.28) 

^"0 = ^"0,0 + ^"0,1 All + t'0,2A22 + o(Aii, A22), 

c = Co + ciAii + C2A22 + o(Aii, A22)- 

Expanding functions "djjj — 1, 2, 3, 4 in the system (5.23), together with substitution 
of assumption (5.28), the second and third equation is divided by An and A22, respec- 
tively; the fourth equation is divided by A11A22, and letting vqq — 0, An, A22 — > 
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we then obtain 



co = 



- Svrwa.o + 327r^a^ = 0, 

[Sn'^itui^o + W2,o) + 327r^(«i + a2)^]Ai2 



-1 _ 

12 ~ 



0. 



(5.29) 



The first three equations in the system (5.9) have a solution 



Co = 0, wo,o = 0, ujifi 



1^2,0 = 47r^Q;2i 



(5.30) 



The fourth equation in the system (5.29) satisfied automatically by using (5.25) and 
(5.30), thus we have 



Using (5.28) and (5.31), we conclude that 

fo = o(Aii, A22) — ^ 0, c = o(Aii,A22) — >0, wi = 47r^a^ + o(Aii, A22) — > 4:7i^al, 

= 471^^2 + o(Aii, A22) — y An'^al, as Aii,A22-^'0, 

and therefore we have (5.26). So the two-periodic wave solution (5.23) tends to the 
two-supersoliton solution (5.6) as An, A22 — )• 0. □ 

6. Conclusion and future work 

Following the procedure described in this paper, we are able to construct quasi- 
periodic wave solutions for other supersymmetric equations also can be dealt with 
by the same way. For instance, 

(1) Supersymmetric KdV equation [311 SH 112] 



Co = Ci = C2 = 0, fo,o = 0, Uifl 



(5.31) 



't + 3 ($2)$)^ + $, 



XXX 
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(2) Supersymmetric MKdV equation [lO lHSl 161] 

ct)t - 30D(0,)S)0 - 3(D0)Vx + <Pxxx = 0, 

(3) Supersymmetric It's equation |65] 

^tFt + QiF,iT)tF)), + ^tF,,, = 0, 

(4) Supersymmetric two-boson equation [66], [67] 

= 2)((D0l)2)+202,x-01,xx, 
4>2,t = 2((I)0l)02)x + 4>2,xx- 

The system (3.10) indicates that constructing multi-periodic wave solutions de- 
pends on the solvability of the system. We consider the number of constraint equa- 
tions and some unknown parameters. Obviously, the number of constraint equations 
of the type (3.10) is 2^. On the other hand we have parameters th, Tij, ai, Ui, (j)o, c, 
whose total number is ^N{N + 1) + 3A^ + 2. Among them, 2N parameters tu, Ui are 
taken to be the given parameters related to the amplitudes and wave numbers of A^- 
periodic waves. Hence, the number of the unknown parameters is |A^(A^+l) + A^+2. 
while |A^(A^+1) parameters Tij implicitly appear in series form, which is general can 
not to be solved explicit. Hence, the number of the explicit unknown parameters 
is only + 2. The number of equations is larger than the unknown parameters 
in the case when N > 2. This fact means that if equation (3.10) is satisfied by 
the unknowns, we have at least A^-periodic wave solutions (A^ < 2). It is still pos- 
sible to construct multi-periodic wave solutions for A^ > 3 by adding the number 
of parameters (for example, using constant solution and integration constant) or 
decreasing the number of equations (for example, using odd and even properties of 
considered equations). In this paper, we consider one-periodic wave solution of the 
equation (1.2), which belongs to the cases when A^ = 1 and A^ = 2 in the Riemann 
theta function (3.1). There are still certain difficulties in the calculation for the case 
A^ > 3. 
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We expect the proposed method to be extended to J\f — 1 supersymmetric 
sine-Gordon equation and M — 1 supersymmetric KP equation, as well as other 
discrete supersymmetric equations (like supersymmetric Toda lattice). For the M = 
2 supersymmetric equations with ordinary variables x, t and Grassmann variables 
61,62, their corresponding superspace is R^'^ — A^xAf. In this case, the Grassmann 
algebra G2(ci, (72) whose dimension is four. The r matrix will be dependent on the 
odd parameters a"i,(j2. In superspace M^'^, the super bilinear forms oi Af = 2 
supersymmetric equations, their quasi-periodic solutions and asymptotic properties 
remain open. We intend to return to these question in some future publications. 
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